Quantum illumination with Gaussian states 
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An optical transmitter irradiates a target region containing a bright thermal-noise bath in which 
a low-reflectivity object might be embedded. The light received from this region is used to decide 
whether the object is present or absent. The performance achieved using a coherent-state transmitter 
is compared with that of a quantum illumination transmitter, i.e., one that employs the signal 
beam obtained from spontaneous parametric downconversion (SPDC). By making the optimum 
joint measurement on the light received from the target region together with the retained SPDC 
idler beam, the quantum illumination system realizes a 6 dB advantage in error probability exponent 
over the optimum reception coherent-state system. This advantage accrues despite there being no 
entanglement between the light collected from the target region and the retained idler beam. 

PACS numbers: 42.50.Dv, 03.67.Hk, 03.67.Mn 
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Entanglement is a fundamental quantum mechani- 
cal resource, with potential applications to Heisenberg- 
limited precision measurements teleportation 0, and 
quantum cryptography Q. Loss and noise, however, 
can quickly destroy this entanglement. Recently, Lloyd 
building on the work of Sacchi 5], has shown that 
"quantum illumination" can reap substantial benefits, 
from the use of entanglement in target detection, de- 
spite the presence of entanglement-destroying loss and 
noise. In Lloyd's quantum illumination paradigm, a pho- 
tonic source creates d-mode maximally-entangled signal 
and ancilla beams each containing a single photon. The 
signal beam irradiates a target region containing a very 
weak thermal-noise bath — with an average of & <C 1 pho- 
tons per mode — in which a low-reflectivity (r; <C 1) ob- 
ject might be embedded. The light received from this 
region — together with the retained ancilla beam — is then 
used to decide whether the object is present or absent. 
Lloyd showed that quantum illumination, with the op- 
timum joint measurement on the received light and the 
ancilla, achieves an effective signal-to-background ratio 
of rjd/b, whereas optimum quantum reception of light 
received in response to transmission of a single unen- 
tangled photon has a much lower (for d 3> 1) signal- 
to-background ratio of rj/b. 

The analysis in was confined, primarily, to the vac- 
uum plus single-photon manifold, wherein at most one 
photon arrives at the receiver during the measurement 
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interval, regardless of whether the object of interest is 
absent or present in the target region. Lloyd did briefly 
describe a microcanonical noise model whose quantum- 
illumination performance — for multi-photon signal and 
noise states with the same energy — again showed the 
factor of d improvement in the signal-to-background ra- 
tio. However the microcanonical model is rather restric- 
tive, in that thermal-noise baths arc in Gaussian states 
with Bose-Einstein distributed photon numbers. In this 
Letter we will remedy that deficiency by providing a 
full Gaussian-state treatment of quantum-illumination 
target detection. Our treatment will employ the ex- 
act quantum statistics for the entangled signal and idler 
beams obtained from continuous-wave (cw) spontaneous 
parametric downconversion (SPDC) in the absence of 
pump depletion [6]. It will also use the standard model 
for the lossy, noisy bosonic channel 0. We will show 
that in a very lossy, very noisy environment, a low- 
brightness quantum illumination system enjoys a sub- 
stantial improvement in effective signal-to-background 
ratio — which can translate into a very large reduction of 
the target-detection error probability — in comparison to 
that achieved by a coherent-state transmitter of the same 
average photon number. As seen in this performance 
advantage accrues to the quantum illumination system 
despite there being no entanglement between the light 
that is received from the target region and the retained 
idler beam. Quantum illumination thus becomes the first 
example of an entanglement-based performance gain, in a 
full bosonic-channel setting, that survives entanglement- 
killing loss and noise. 

Consider an entangled signal and idler mode pair ob- 
tained from cw SPDC. This mode pair — with annihila- 
tion operators ds and di — is in the entangled state with 
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number-ket representation 
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where Ns is the average photon number per mode. In the 
quadrature representation, \'ip)si is a zero- mean Gaussian 
state whose Wigner-distribution covariance matrix is 
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where S = 2Ns + 1 and Cq = 2^Ns[Ns + 1). Viewed 
in this way, it is easily seen that \'4>)si has maximally- 
entangled quadratures, because the magnitudes of non- 
zero ofF-diagonal terms, Cq, in A.si equal the maximum 
value allowed by quantum mechanics given the diago- 
nal elements of that matrix. Indeed, the upper limit on 
the magnitudes of these off-diagonal terms for a classi- 
cal state, i.e., one with a proper P representation so that 
the signal and idler modes are unentangled, is Cc = 2A^s. 
Thus, for a low-brightness source, for which the average 
number of photons per mode is very low {Ns ^1), there 
is a very strong nonclassical signature in A5/ in that 
Cq » Cc prevails. 

Suppose that we transmit the signal mode toward a 
spatial region that may or may not contain a weakly- 
reflecting target but, in either case, contains a bright 
thermal-noise bath. Also suppose that the transmitter 
retains the idler mode, for a subsequent joint measure- 
ment to be made with the return from that target region. 
Under hypothesis Hq (target absent), the annihilation 
operator for the return from the target region will be 
Sfl = SiB^ where as is in a thermal state with average 
photon number Nb ^ 1. Under hypothesis Hi (target 
present), the return- mode's annihilation operator will be 
0.R = v^^s + Vl ~ kob, where k and as is now in a 
thermal state with average photon number Nb/{1 — n)- 
Physically, this represents a very lossy (k ^ 1) return 
from the target when it is present, combined with a very 
strong background contribution {Nb >• 1) that, at the 
receiver, is independent of target absence or presence Q. 

Under both hypotheses the joint state of the qr and di 
modes is a zero-mean mixed Gaussian state, with condi- 
tional Wigner-distribution covariance matrices given by 
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under ffi, where B = 2Nb + 1 and A = 2kNs + B. 
Neither of these conditional states is entangled. This is 
clearly so when the target is absent, because the bath 
mode Ob is independent of the idler mode a/. When the 
target is present the lack of entanglement when Nb > k 
follows from the magnitude of the off-diagonal elements, 
y/nCq, in A^] falling below the classical-state limit, 
2 \/Ns{kNs + Nb), set by the diagonal elements in that 
matrix. In other words, the conditional density opera- 
tors, /5^] and p^j, for the return and idler modes under 
hypotheses Hq and Hi both have proper P representa- 
tions, and hence can be regarded as classically-random 
mixtures of coherent states. 

When the two hypotheses are equally likely, the min- 
imum error probability decision rule for the quantum il- 
lumination receiver is as follows Q . Measure p^] — p^] 
and declare the target to be present if a non-negative 
outcome results, and otherwise declare the target to be 
absent. Moreover, with {7^^^} denoting the non-negative 
eigenvalues of p^] — p^] , we have that the error probabil- 
ity of this optimum quantum receiver for the quantum- 
illumination transmitter is 
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Unfortunately, finding the eigenvalues of p^] — p^] is 
a daunting task, so, as is often done in communication 
theory, we shall establish bounds on this error probabil- 
ity. Before doing so, however, we pause to introduce two 
classical-state comparison cases for the quantum illumi- 
nation system. The first uses a coherent-state transmit- 
ter with average photon number Ns- The second em- 
ploys signal and idler modes that are in a joint classical 
state with only the signal being transmitted toward the 
target region, while the idler is retained for possible use 
in conjunction with return mode from the target region. 
Optimum quantum reception will be assumed for both 
of these cases, so that they too employ difference of con- 
ditional density operator measurements and have error 
probabilities given by in terms of the non-negative 
eigenvalues of their measurement operators. 

Based on [4] , which showed that quantum illumination 
is advantageous at low signal-to-background ratios, and 
the fact that the off-diagonal elements in A5/ asymp- 
totically approach classical-state behavior for Ns ^ 1, 
we have chosen to study quantum illumination in the 
regime wherein Ns <^ 1 and Nb ^ 1. Only a moment's 
thought is needed to conclude that the error probability 
in this operating regime will be very close to 1/2 for our 
quantum illumination system, as well as for the coherent- 
state and classical joint-state comparison cases. So, to 
get to acceptably low error probabilities with the quan- 
tum illumination and coherent-state systems we will use 
Af 3> 1 identical transmissions of the types described 
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above, in conjunction with optimum joint quantum mea- 
surements. Thus, the joint signal-idler state for the quan- 
tum illumination transmitter will be pgj = <8)m=iPSm-fm ' 
where the modal signal-idler states are zero-mean and 
jointly Gaussian with the Wigner distribution covariance 
matrix from and the coherent-state transmitter will 



emit 'S>m=l\^/Ns)m, where {|-)m} are the modal coher- 
ent states. The description of our general classical-state 
transmitter will be given later. 

At this point, the quantum Chcrnoff bound [lo| comes 
to our rescue. For optimum quantum discrimination be- 
tween a pair of equally likely Af-mode conditional density 
operators, p'*^-* — ®m=iPrn for A: = 0, 1, with identical 
modal states under each hypothesis, viz., p'm = p^'^^ for 
1 < m < M and k = 0, 1, the quantum Chernoff bound 
places the following limit on the error probability: 



Pr(e) < le-''^ = - 



min trf (p^^^Yip'^^H^- 
0<s<l V 



n M 
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This bound is exponentially tight, i.e., the error probabil- 
ity exponent, — ln[Pr(e)]/M, converges to £ as M oo. 
Its potentially weaker s = 1/2 version, known as the 
Bhattacharyya bound. 
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will also be of interest in what follows, because it is more 
amenable to obtaining analytic results and because it is 
related to the lower bound on the error probability 

Pr(e) > ^ (l - ^1 - [tr((p(0))i/2(p(i))i/2)]2^^ 

which, in general, is not exponentially tight. 

Our quantum illumination and coherent-state trans- 
mitters both lead to modal density operators, under 
each hypothesis, that are Gaussian states. Hence we 
can use the results of [lH to evaluate Chernoff or Bhat- 
tacharyya bounds. To do so we need three things: the 
mean values of the relevant modes — for the coherent- 
state transmitter, and clr and d/ for the quantum illu- 
mination transmitter — under each hypothesis; the condi- 
tional Wigner-distribution covariance matrices for these 
relevant modes; and the symplectic diagonalizations of 
those conditional Wigner-distribution covariance matri- 
ces. The symplectic diagonalization of a 2K x 2K di- 
mensional covariance matrix A consists of a 2K x 2K 
dimensional symplectic matrix S and a symplectic spec- 
trum {vk : 1 < k < K} that satisfy 
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where diag(-, •,...,•) denotes a diagonal matrix with the 
given diagonal elements. 

For the coherent-state transmitter we have that 
tr(p(o)a^i) = 0, tr(p(i)afl) = 03Vs, A(°) = A^^) = 
diag(B/4,E/4). It follows that S is the 2-D identity 
matrix and vi = B/A, so that the results in [ll[ lead to 
the quantum Chcrnoff bound (which turns out to be the 
Bhattacharyya bound) 

Pr(e)cs < e-^'^'^^i^^^^-^y /2 (11) 
_ g^M«Ars/4iVB/2, when A^B > 1. (12) 
Because (|11|) is also the Bhattacharyya bound, we have 
Pr(e)cs > 

i ^1 _ ^/l - e-'^MnNsiVN^-Vm)'''^ (13) 

when A^B > 1 and MkNs/2Nb > 1. 

For the quantum-illumination transmitter, the 4-D 
identity matrix is the symplectic matrix needed for the 
diagonalization of A^] , and 
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is the symplectic matrix needed to diagonalize A^]. 




Here, X± = diag(2;±, ±a;±) with 
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The associated symplectic spectra are vi = B/A,V2 
5/4 under Hq, and 



{-lf{S~A) + J{A + S) 



/8, (17) 



A = Sdia,g{yi,yi,V2,^'2,---,VK,i'K)S^,{lO) 



for k — 1,2, under Hi. These diagonalizations can be 
employed to derive an analytic expression for the Bhat- 
tacharyya bound on the error probability achieved with 
quantum illumination. Unfortunately, that expression is 
far too long to exhibit here. In Fig. 1 we compare the 
coherent-state system's Chernoff bound from pT|) and 
its lower bound from (jl3p with the quantum illumina- 
tion system's Bhattacharyya bound when k — 0.01, Ns = 
0.01, and Nb — 20. We see that the quantum illumina- 
tion system's error probability upper bound — at a given 
M value — can be orders of magnitude lower than the er- 
ror probability lower bound for coherent-state light. 

To show that the advantage afforded by quantum illu- 
mination extends well beyond the specific example cho- 
sen for Fig. 1, we have used an algebraic computation 
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FIG. 1: Upper bounds (solid curves) on the target-detection 
error probabilities for coherent-state (Chernoff bound) and 
quantum illumination (Bhattacharyya bound) transmitters 
with M transmitted modes each with Ns ~ 0.01 photons 
on average when k — 0.01 and Nb = 20. Also shown is the 
lower bound (dashed curve) for the coherent-state case, which, 
see below, also applies to all classical-state transmitters with 

program to obtain the following approximate form for 
the quantum illumination transmitter's Bhattacharyya 
bound when k < 1, TV^ < 1 and Nb > 1: 

Pr(e)Qi < e-*^"^«/^«/2. (18) 

Comparing this low-brightness quantum illumination re- 
sult to the Chernoff bound for the coherent-state trans- 
mitter in the same lossy, noisy environment reveals that 
the use of entanglement at the transmitter improves the 
error probability exponent by a factor of 4, i.e., by 6dB. 
One might wonder whether there is some other classical- 
state source that can match the performance of the quan- 
tum illumination transmitter in this regime. To show 
that such is not the case, we now develop the perfect- 
measurement bound on the performance of any classical- 
state transmitter. Consider the general classical-state 
(signal and idler) transmitter whose output state is 

Psi = J Jd^asd^aj P{as,ai)\as)s\cti)is{as\i{ai\, 

(19) 

where {|Q:j)j} for j = S,I are the M-mode coherent 
states of the signal and idler and ^(0:5,0:/) is a prob- 
ability density function over {aj : j — S, I } satisfy- 
ing J2m=ii^S ^-Sm) — MNs- Minimum error probabil- 
ity reception when the signal beam irradiates the target 
region and a joint measurement is made of the return 
from the target region and the retained idler is bounded 
from below by the error probability achieved, using this 
transmitter, when 0:5 is known to the receiver. A noisy 
estimate of can be obtained from heterodyne detec- 
tion of cy.1, hence our lower bound is correctly termed a 
"perfect-measurement" bound. More importantly, given 



knowledge of as, the classical-state system reverts to be- 
ing a coherent-state transmitter, so that (fT3|) . convexity, 
and the average photon-number constraint imply that 
Pr(e)ciass obeys the same lower bound as Pr(e)csj show- 
ing that when Nb ^ 1 all classical-state transmitters 
have error probability exponents that are at least 3 dB 
inferior to that of quantum illumination with the same 

Em=i(aL"-5™) value. 

Given that Nb = k suffices to destroy the entangle- 
ment between the as^ and an^ modes under hypothesis 
Hi, it behooves us to find a physical explanation for the 
performance gain provided by quantum illumination. In 
our view the explanation is as follows. When Ns <C 1, 
the entangled as^ and aj^ modes have a cross correla- 
tion that greatly exceeds what is permitted for a classical 
state. Although the corresponding cross correlation for 
dij^ and d/^ when the target is present is within clas- 
sical bounds, there is no classical input state for ds„ 
and d/^^ of the same average photon number that can 
produce a close approximation to this output state. Of 
course, a single- mode transmission with Ns ^ 1, k 1, 
and Nb 3> 1 cannot yield a low error-probability de- 
cision. However, the joint measurement over Af 3> 1 
such modes can. Moreover, no special joint-state perfor- 
mance enhancement is realized with the coherent-state 
transmitter, because the product state (8)^=1 1 V^s)™, is 
merely a coherent state \^/MNs) of a different mode. In 
summary, we must be careful not to dismiss the value 
of using entangled resources just because the application 
scenario is lossy and noisy. 
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